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ABSTRACT 


In a continuous search for the energy-efficient electronic switches, a great attention is focused on tunnel field-effect transistors 
(TFETs) demonstrating an abrupt dependence of the source-drain current on the gate voltage. Among all TFETs, those based 
on one-dimensional (1D) semiconductors exhibit the steepest current switching due to the singular density of states near the 
band edges, though the current in 1D structures is pretty low. In this paper, we propose a TFET based on 2D graphene bilayer 
which demonstrates a record steep subthreshold slope enabled by van Hove singularities in the density of states near the 
edges of conduction and valence bands. Our simulations show the accessibility of 3.5 x 10"' ON/OFF current ratio with 150 
mV gate voltage swing, and a maximum subthreshold slope of (20 /xV/dec)^' just above the threshold. The high ON-state 
current of 0.8 mAZ/xm is enabled by a narrow (~ 0.3 eV) extrinsic band gap, while the smallness of the leakage current is due 
to an all-electrical doping of the source and drain contacts which suppresses the band-tailing and trap-assisted tunneling. 


Introduction 

The design of field-effect transistors (FETs) operating at sub-0.5 V supply voltage is one of the major challenges for nano¬ 
electronics paving the way to resolve the problem of power dissipation in large integrated circuits. Tunnel FETs (TFETs) 
with interband tunneling are among the principal candidates to meet this demand.'’^ The low-voltage switching in TFETs 
occurs due to a sharp dependence of the tunnel current on the conduction-valence band overlap in a gate-controlled p-i or 
p-n junction;^ once there is no band overlap, there is no tunnel current. This fact, alongside with the smallness of thermionic 
leakage current and efficient modulation of the barrier transparency by the gate voltage, results in the subthreshold slope of 
the TFET characteristics surpassing the thermionic limit"'"" of (60 mV/dec)"*. 

It is intuitive that the abrupt variations of electron and hole densities of states (DoS) near the band edges further enhance 
the switching efficiency of TFETs. ^ For a x/-dimensional TFET channel, the DoS scales with energy E above the band edge as 
p;{d- 2 )l 2 ^ while the dependence of current J on the gate voltage Vq above the threshold Vt is J [Vq — " Apart 

from the density-of-states-enhanced switching, the TFETs with low-dimensional channels demonstrate an improved electro¬ 
static control of the band alignment by the gate voltage."’ Theoretically, the effects of DoS on current switching steepness 
are most pronounced in vertical TFETs based on the two-dimensional crystals'"*’'^ and electron-hole bilayers in quantum- 
conhned structures.'^’In such TFETs, the joint density of states is nonzero just at one certain value of gate voltage'" 
- which could lead to the abrupt-most current switching ever. In practice, the density-of-states effects on the subthreshold 
steepness are largely smeared. The reason for the smearing in vertical TFETs based on van der Waals heterostructures is the 
rotational misalignment of 2D layers in momentum space.'® In common semiconductor structures it is believed that charged 
defects and dopants lead to pronounced band-tailing and emerging trap-assisted and band-tail tunneling leakage currents. 

In this paper, we theoretically demonstrate that graphene bilayer (GBL) represents an ideal platform for the low-voltage 
tunnel switches. A peculiar ’mexicanhat’ band structure of GBL formed under transverse electric field^"’^"' leads to a van Hove 
singularity in the DoS right at the bottom of the band edges, as shown in Fig. 1 A. It was experimentally proved that this singu¬ 
larity leads to an enhanced infrared absorption.^^ Another manifestations of this singularity predicted theoretically include a 
step-like change of the GBL FET channel conductance^® and multi-peak structure of the vertical GBL TFET characteristics.^^ 
In this paper we show that such singularity results in a steep, linear dependence of the GBL TFET current on the gate voltage 
above the threshold, which was attributed previously just to the TFETs based on one-dimensional materials.®’"’^"""" 

The advantage of graphene bilayer TFETs over those based on 2d materials with parabolic bands in terms of switching 
steepness can be illustrated by Fig. 2 B. As the conduction and valence bands in a GBL tunnel junction overlap, the electrons 




Figure 1 . (A) Electron spectrum E[p) in graphene bilayer under transverse electric field and the energy dependence of its 
DoS. The ’’Mexican hat” feature in the dispersion law leads to the square-root singularities in the DoS near the band edges. 
Panel (B) highlights with red the electron states involved in the interband tunneling at small band overlap in graphene bilayer 
(left) and in a semiconductor with parabolic bands (right). The phase space for tunneling in graphene bilayer represents a 
ring, while in a parabolic-band semiconductor it is a point. Dashed lines indicate the tunneling transitions, red lines indicate 
the trajectories of the tunneling electrons in the valence band. 


capable of tunneling are located on a ring in the momentum space. In contrast, the electrons capable of tunneling between 
simple parabolic bands are located in a small vicinity of an extreme point of the dispersion. To fully exploit the density-of- 
states effect for tunneling and get rid of the dopant-induced band edge smearing, we introduce the TFET based on electrical 
doping by auxiliary gates. Apart from avoiding the parasitic tunnel currents, this adds the possibility to electrically reconfigure 
the device between n- and p-types. Under optimal gate biasing conditions, the proposed TEET demonstrates the current 
switching over more than 4 orders of magnitude with 150 mV gate voltage swing only. At the same time, the ON-state current 
density as large as 0.8 mA//rm is accessible due to the low extrinsic band gap of GBL (~ 0.3 eV) and large DoS far above the 
band edges. 

Results 

Device structure 

The advantages of graphene bilayer for the steep current switching can be fully realized in the structure of the TFET shown 
schematically in Fig. 2 A. A heavily doped silicon substrate acts as a bottom gate used to create the transverse electric field 
and thus open and manipulate the band gap in GBL.^* The oxidation of the substrate results in formation of Si02 layer playing 
the role of back gate oxide and substrate for graphene bilayer. Alternatively, the Si02 layer can be replaced with hexagonal 
boron nitride possessing a small (~ 10**^ cm^^) density of residual charged impurities.A nanometre - thin layer of high-ff 
dielectric (e.g., zirconium oxide) covers the graphene channel, and the top metal gates are formed above. The side gates near 
the source and drain contacts induce large densities of holes and electrons, respectively, which also leads to the formation of 
an abrupt tunnel junction and energy barriers (see Fig. 2 B) for the thermally activated electrons and holes contributing to the 
OFF-state leakage current. 

The operation of a normally open TFET switched off by a negative top gate voltage is illustrated in the band diagram. 
Fig. 2 B. Application of positive voltage to the bottom gate, Vb > 0, induces the band gap and provides an excess electron 
density in bilayer. The /?+ doping of source emerges upon application of negative voltage f/j < 0 to the source doping gate. An 
additional increase in the barrier height for the holes injected from the drain is achieved by applying positive voltage Ud > 0 
to the drain doping gate. It is instructive that application of high voltage to the doping gates does not result in increased 
power consumption as this voltage is not changed during the device operation. At zero top gate voltage, the valence band in 
p+-source overlaps with the conduction band in the n-type channel, which corresponds to the ON state (red band profiles 
in Fig. 2 B). Upon application of negative voltage to the top gate, the transistor is switched off (dashed blue band prohles in 
Fig. 2 B). 

The optimization of the device dimensions aiming at the increase in the ON-state and reduction in the OFF-state currents 
is quite straightforward: both the effective thickness of the gate dielectric and the distance between the source doping gate and 
the control gate should be small. These distances are limited just by the possible gate leakage current (see below), we choose 
them to he dt =2 nm and dg = 5 nm. The doping gate at the drain is used just to induce high barrier for thermally activated 
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Figure 2. (A) Layout of the proposed graphene bilayer TFET with electrically defined source and drain regions (B) Band 
diagram of graphene bilayer TFET for the optimal biasing conditions; Vb > 0,11$ < 0, Ud > 0. At zero top gate bias, Vq = 0, 
the TFET is switched on, while at Vg < 0 it is switched off 


holes; the distance between this gate and the control gate should be large to reduce the transparency of tunnel junction at the 
drain and get rid of ambipolar leakage. 

Model of the interband tunneling enhanced by van Hove singularities 

Our modeling of graphene bilayer TFET relies on a self-consistent determination of the carrier density and band structure^^ 
followed by the calculation of tunnel current under assumption of ballistic transport (see supplementary material, sections I- 
111). However, the principal dependence of the tunnel current on the gate voltage can be derived in a very simple fashion. The 
current is proportional to the number of electrons capable of tunneling between the valence band of source and the conduction 
band of channel. Once these band overlap by dE in the energy scale, the electrons available for tunneling in graphene bilayer 
occupy a ring in the momentum space (Fig. 1 B, left panel). Their number is proportional to pi-amdE, where pmin is the 
momentum corresponding to the bottom of the ’Mexican hat’. One thus concludes that the tunnel current is a linear function 
of the band overlap which, in turn, is a linear function of the gate voltage. This contrasts with the 2d materials having parabolic 
bands where the number of electrons available for tunneling is proportional to \/EdE (Fig. 1 B, right panel). As a result, the 
current in TFETs based on these materials is proportional to the gate voltage raised to the power 3/2. 

A rigorous expression for the tunnel current density involves an integral of the single-particle velocity V|| = dE/dp^\ 
timed by the barrier transparency &{p±,E) and the difference of occupation functions in the valence and conduction bands 
/y(£’) — fc{E) over the momentum space d^p = 2dp±^dp^\ 

Ey Pm‘dx{E] 

Ji=^JdE[fAE)-ME)] J 2ufpxx2^(£,px)- (D 

Ec 0 

Here, gsgv = 4 is the spin-valley degeneracy factor in graphene, pmaxiE) is the maximum transverse momentum of electron at 
a given energy E, pmax(F) = inin{pc(F),Pv(F)}, where pc{E) and Pv{E) are the inverse functions to the electron dispersion 
in the conduction and valence bands. The limits of integration over energy are the conduction band edge in the channel, Ec, 
and the valence band edge in the source, Ey. The factor of two before the quasi-classical barrier transparency comes from the 
presence of two turning points with zero group velocity in the GBL dispersion at which an electron attempts to tunnel. 

The effect of ’Mexican hat’ on the current switching steepness can be traced analytically from Eq. 1 by assuming that 
the conduction band states are empty, valence band states are occupied, and the barrier transparency ^{E,p±) ^ S’q weakly 
depends on the energy and transverse momentum. At small band overlap, the momenta of the tunneling electrons in graphene 
bilayer are close to pmin (Fig- IB, left panel), which results in 

( 2 ) 

This linear dependence is in agreement with the above qualitative considerations. Previously, such a dependence of the 
tunnel current on the band overlap was attributed just to the ID semiconductor structures which proved to be among the best 
candidates for the TFETs. 

An additional increase in the graphene bilayer TFET subthreshold steepness occurs due to the dependence of the trans¬ 
parency 3l{E,pi_) on the junction field and, hence, gate voltage. The transparency is evaluated by integrating the imaginary 
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Figure 3. Calculated room-temperature gate transfer (left) and current-voltage (right) characteristics of graphene bilayer 
TFET at hxed bias voltages at auxiliary gates: Vg = 3.3 V, Us = —0.6 V, Ud = 0.25 V. Top gate dielectric is 2 nm Zr 02 , 

K = 25, back gate dielectric is 10 nm Si 02 , spacing between the gates dg = 5 nm. Inset: gate transfer characteristic in the log 
scale 


part of the electron momentum inside the band gap, which results in (see Supporting information, section III) 


^(£',px)~exp 


-|lmp||(£ = 0) 


(3) 


where Imp|| {E = 0) is the imaginary part of electron momentum evaluated at the midgap, and I is the length of the classically 
forbidden region (tunneling path length). The latter is given hy I = A/ eF for px < Pmm, I = 2E {pi_)/eF for px > Pmm^ 
where F is the electric held at the junction found from the solution of Poisson’s equation, and A is the band gap in the GBL. 
The thermionic leakage currents were evaluated with equations similar to (1) by setting the unity transmission probability and 
constraining the energy integral to the particles with the energies above the barrier. 


Characteristics of the graphene biiayer TFET 

The calculated room-temperature /(Vcj-characteristics of graphene bilayer TFET at different drain bias Vd are shown in Fig. 3 
A. The current density just above the threshold voltage is a linear function of Vq, in agreement with the simple density-of- 
states arguments and Eq. (2). With increasing the top gate voltage, the J{Vg) - curve becomes superlinear, which is attributed 
to the exponential sensitivity of the tunnel barrier transparency to the junction field. The subthreshold slope at Vq = Vth 
reaches (20 pV/dec)^^ and is limited by the small thermionic current Jff, ~ 1.3 x 10^^ mA/pm and the gate leakage current 
~ 1.0 X 10^^ mA/pm. However, it is not the highest subthreshold slope that determines the power efficiency of the TFET 
but rather the supply voltage V 5 required to switch the transistor between the ON- and OFF-states. Considering the current at 
Vg = 0 V as the ON-state current (Jqn = 0.8 mA/pm in Fig. 3 B at Vo = 0.15 V), and the leakage current as the OFF-state 
current (Jon/Joff = 3.5 x 10^), we have obtained V 5 = 150 mV. In a conventional MOSFET, the gate voltage swing V 5 > 285 
mV is required to achieve the same current switching ratio. 

The drain characteristics of graphene bilayer TFET shown in Fig. 3 B demonstrate a pronounced current saturation typi¬ 
cally absent in single graphene layer FETs. This saturation is due to the limited energy range in which the tunneling injection 
is possible. The presence of saturation is important for the logic inverters which guarantees the clear discrimination of the 
zero and unity output signals. 

The average subthreshold slope of our TFET over 4.5 decades of current is 33 (mV/dec)^^ With this characteristic, it 
outperforms all sub-thermal tunnel switches^ based on silicon,"^ germanium,® III-V hetero junctions,® and carbon nanotubes^ 
reported to date. Only recently a vertical TFET based on MoS 2 /germanium junction with a similar value of the average 
subthrehold slope was demonstrated,®^ however, its ON-state current density of 1 pA/pm leaves much to be desired. 

The aggregate quality of the TFET, accounting for both average subthreshold slope and current density, can be char¬ 
acterized by an /go-figure of merit®® which is the current density at the point where the subthreshold slope equals (60 
mV/dec)^^ While the best /go reported to date equals 6 nA/pm (InAs nanowire/Si heterojunction TFET®), in our TFET 
structure /go =150 pA/pm. 
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Figure 4. Comparison of the gate transfer characteristics of GBL TFET and a TFET based on an equivalent 2D 
semiconductor with the same barrier transparency but with different (parabolic) band structure. Numerical values of the 
effective masses are taken for bulk InAs. The insets show the band diagrams overlaid with electron-hole spectra and the 
energy dependence of DoS 


The unique characteristics of GBL TFET surpassing the existing TFETs are enabled by the three factors. First of all, it is 
the small extrinsic band gap (for doping gate voltages used in Fig. 3, A Ri 0.3 eV) that guarantees elevated tunneling probability 
(^ ^0.1) and large current density. Most two-dimensional transition metal dichalcogenides have large intrinsic band gaps 
(1.9 eV for M 0 S 2 , 1.3 eV for WS 2 , etc.), while in the 2D structures based on III-V materials being narrow-gap in the bulk, the 
gap value increases significantly due to the quantum confinement.^® Secondly, the singular DoS near the band edges allows an 
abrupt switching of tunnel current. Even if there existed a parabolic-band 2D material with the same band gap and the same 
barrier transparency in the TFET structure, its current density would be given by (see supplemental material, section IV) 




par ' 



Imctriv 
nic -f m, 




(4) 


where rric and m,, are the conduction and valence band effective masses and, similar to the derivation of Eq. (2), we have 
assumed the barrier transparency to be energy- and momentum independent. The numerical comparison of current density 
in graphene bilayer and its equivalent parabolic band counterpart is presented in Fig. 4 for the effective mass values typical 
for narrow-gap III-V semiconductors (me = 0.024mo, = 0.026mo for InAs). At 150 mV gate voltage above the threshold, 

the current density in graphene bilayer exceeds 15 times that in a parabolic-band material. The factor of two is due to the 
valley degeneracy absent in III-V’s, another factor of two is due to the tunneling at two turning points of the ’Mexican hat’ 
dispersion, and the remainder of 3.5 is due to the finiteness of electron momentum at the edge of the ’Mexican hat’. 


Gate leakage and band tailing: the insulator selection rules 

The steep switching of the tunnel current by the gate voltage can be masked by the leakage to the gates, band-tail and trap- 
assisted tunneling.The latter factors might have masked the onset of the interband current in the recent measurements of 
graphene bilayer tunnel junctions.A careful selection of the gate dielectrics providing high interface quality is required 
to minimize these effects. 

The main reason for the band tailing comes from the fluctuations of electric potential produced by the random charged 
defects or dopants.^ This effect is most pronounced in the TFETs with source and drain intentionally doped chemically. In 
the TFETs with electrically doped contacts, only residual charged impurities inevitably present in the substrate contribute 
to the band tailing. To provide a quantitative view on the band tailing in graphene bilayer on different substrates, we have 
evaluated the quasi-classical DoS p{E) in the presence of fluctuating potential by integrating the singular ’bare’ DoS po over 
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Figure 5. Calculated energy dependencies of the DoS in the conduction band of graphene bilayer at different densities of 
charged impurities (corresponding to the substrates of different quality). The electron density is held hxed at 4 x 10*^ cm^^, 
the nominal energy gap is 0.3 eV 


the probabilities of voltage fluctuations^ 


“ Tim / (- 5 m ) 


(5) 


where (V^) is the root-mean-square amplitude of the voltage fluctuations proportional to the impurity density n,. The cal¬ 
culated energy dependencies of the ’smeared’ DoS are shown in Fig. 5. For the parameters of chemical doping used in the 
pioneering proposal of the GBL TFET, n, = 4 x 10*^ cm^^, the conduction and valence bands almost merge together, which 
would result in a poor OFF-state, nothing to say about high switching steepness. A slight peak in the DoS near the band bot¬ 
tom becomes noticeable already at impurity density of 5 x 10^^ cm^^ which corresponds to the low-quality graphene on Si02 
substrates. In graphene samples on a high-quality Si02,^*' the smearing of the band edge is order of 10 meV. The ultimate band 
abruptness of ^ 5 meV can be achieved in graphene samples encapsulated in boron nitride, providing the residual impurity 
density of ^ 5 x 10**' cm^^.^^ At this limit, the fluctuation-induced smearing of the bands becomes negligible, and the behav¬ 
ior of the DoS near the bottom of the ’’Mexican hat” is governed by the trigonal warping distortions of electron spectrum due 
to the next-nearest neighbor interactions.^^ Using the exact spectrum of GBL with trigonal warping, we estimate the energy 
scale where the trigonal warping is relevant as 5e ~ 20 meV. Already for relatively small gate voltages, Vo — Vth > de/e, 
these corrections are irrelevant and the linearity of the 7(VG)-characteristic holds. 

The gate leakage may also limit the minimum achievable OFF-state current, while at the same time small effective gate 
oxide thickness is required to efficiently control the band structure in the channel by the gate voltage. Among the common 
high-fc materials, zirconium oxide (k ^ 25) looks as an optimal solution for the GBL TEFT due to the large band offset with 
respect to graphene (Ut = 2.9 eV^*) and elevated tunneling mass m, « 0.3mo.^^ We have evaluated the leakage current from 
graphene with electron (hole) density of into the metal gate to be (see Supporting information, section VI) 


Jg = SeLgii, 


Ub kplio 


f, l + {kFh,g)^ 


% 


g> 


(6) 


where /loc = Ti/y/TmtUb is the electron localization length in the direction perpendicular to the graphene bilayer, kp is the 
Fermi wave vector in the metal, S>g is the transparency of the barrier separating GBL and the gate, and Lq is the gate length. 
Under the biasing conditions of Fig. 3, the gate leakage current is estimated to be 7, = 1.0 x 10^^ mA//zm which is below the 
thermionic leakage level (in this estimate, we have taken Lg = 20 nm and kp =2 A^*). 
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Discussion 


We have proposed and substantiated the operation of a graphene bilayer TFET exploiting the van Hove singularities in the 
density of states near the band edges. The presence of these singularities leads to the increased steepness of the gate charac¬ 
teristics and to the high ON-state current as well. The subthreshold slope of JiVc) curve in the proposed FET reaches the 
maximum of (20 /iV/dec)^^ while only 150 mV gate voltage swing is required to change the current density from Jqn ~ 1 
mA//rm down to Jqff ~ 2 x 10^^ mA//rm. As a matter of fact, the effects of singular DoS on the interband tunneling are 
possible just in the TFET structures with an all-electrical doping, where the effects of band tailing and trap-assisted tunneling 
are minimized. 

Such steep switching in the lateral TFETs based on 2d materials is possible if only the van Hove singularities are present 
both at the top of the valence and the bottom of the conduction band. This property is unique to the graphene bilayer and is 
absent in other 2d materials (e.g., those based on III-V compounds), where a ’Mexican hat’ structure is formed in one of the 
bands due to the spin-orbit coupling.We can thus conclude that graphene bilayer is an only two-dimensional material where 
the switching of interband tunnel current is as steep as in one-dimensional semiconductors, whereas the large on-state current 
is inherited from the single layer graphene. 

Methods 

The modeling of GBL TFET is based on the self-consistent determination of carrier density and band structure under fixed 
gate voltages^^ followed by the calculation of tunnel current with Eq. (1). The necessity for self-consistent calculation is 
dictated by the dependence of the energy gap on the electric field between graphene layers comprising the GBL; the field, 
in turn, depends on the induced carrier density which is sensitive to the band structure. The distribution of electric field 
at the tunneling junction required for the evaluation of the barrier transparency is calculated with the conformal mapping 
technique. The numerical model is described in detail in Supporting information, sections I-II. In section III of the Supporting 
information, an approximate analytic model of GBL TFET is presented. 

The effect of charged impurities present in the substrate on the singular density of states in graphene bilayer is evaluated 
with Kane’s quasi-classical model of band tails. ^ The revision of the model for the two-dimensional GBL is presented in 
Supporting information, section V. 

The gate leakage current is estimated with a quantum-mechanical model of graphene bilayer tunnel coupled to the contin¬ 
uum of delocalized states in the metal gate. The model is presented in Supporting information, section VI. 
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Supporting information to ’’Abrupt current switching 
in graphene biiayer tunnei transistors enabied by 
van Hove singuiarities” 


ABSTRACT 

In a continuous search for the energy-efficient eiectronic switches, a great attention is focused on tunnei fieid-effect transistors 
(TFETs) demonstrating an abrupt dependence of the source-drain current on the gate voitage. Among ali TFETs, those based 
on one-dimensionai (1D) semiconductors exhibit the steepest current switching due to the singular density of states near the 
band edges, though the current in 1D structures is pretty iow. In this paper, we propose a TFET based on 2D graphene biiayer 
which demonstrates a record steep subthreshoid siope enabled by van Hove singularities in the density of states near the 
edges of conduction and vaience bands. Our simuiations show the accessibiiity of 3.5 x 10"^ ON/OFF current ratio with 150 
mV gate voltage swing, and a maximum subthreshoid siope of (20 /xV/dec)^' just above the threshoid. The high ON-state 
current of 0.8 mA/^m is enabied by a narrow (~ 0.3 eV) extrinsic band gap, whiie the smaliness of the leakage current is due 
to an aii-eiectricai doping of the source and drain contacts which suppresses the band-taiiing and trap-assisted tunneiing. 


1 Self-consistent evaluation of band gap and carrier density in gated biiayer 

The intrinsic graphene biiayer (GBL) is a gapless semiconductor with symmetric electron-hole dispersion being parabolic 
near the band edges. However, the application of transverse electric field F± leads to the potential energy difference A = eF±d 
between the layers comprising GBL {d — 3.35 A is the interlayer separation in GBL) and opening of the band gap. As 
the transverse electric field depends on the carrier density, which, in turn, depends on the band structure, a self-consistent 
procedure for the determination of those quantities in gated GBL is required. 

In this section, we describe the procedure used to evaluate the carrier density, bandgap, and electric potential in GBL at 
given potentials of top and bottom gates (V) and Vt). The tight-binding Hamiltonian of graphene biiayer 

/ I v{px-ipy) 7l ® \ 

v{px + ipy) I 0 0 

7i 0 -f v{px + ipy) 

Vo 0 v{px-ipy) -f / 

yields the following energy bands depicted in Fig. 6 :* 
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with the corresponding four-component eigen wave functions (non-normalized): 
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The components of wave function determine the amplitudes of finding an electron on the atoms of type ”A” and ”B” of the 
unit cell on the top and bottom layers in the f-th band. In Eqs. (7) - (10), p is the electron quasimomentum, yi ~ 0.4 eV is the 

















Figure 6. Band structure of graphene bilayer in 

^ is the bandgap, and /?„„■„ = 

+7r 




tight binding approximation. A is the interlayer asymmetry parameter, 
is the electron momentum at the band edges 


interlayer hopping integral, and v Ri 10® m/s is the Fermi velocity in graphene. Given the eigenfunctions, the probabilities of 
finding an electron of the /-th energy band on the top and the bottom graphene layer can be calculated as follows: 


v'i (P) = 


'(P) = 


‘p;""' (p) + 'i’r" (p) 


I'I'Kp) 


Xy^bottom _j_ xy^bottom 


^bottom I 


I'I',(P) 


( 11 ) 

( 12 ) 


here 11... 11 stands for the norm of the vector. 

The net charge densities at top and bottom graphene layers are obtained by summing the probabilities (11) timed by the 
distribution functions over the bands and momenta: 


Pt = eg< Y. ^5(p)[i-/(£i(p):M)]- L ^!(p)/(£/(p),m) 

Ip.'= 3,4 P.j=1.2 


Pb = eg< Y ^?(p)[ 1 -/(^'(p),M)]- Y 

[p,r=3,4 P,!=1,2 J 


(13) 

(14) 


E—ji . 

Here, g = 4 is the spin-valley degeneracy factor, f{E,ji) = [e~'^ + 1] ^ is the Fermi-Dirac distribution function, ji is the 
chemical potential reckoned from the midgap, k is the Boltzmann constant, and T is the absolute temperature. 

We further consider graphene bilayer placed between the top and bottom gates separated by gate dielectrics of thicknesses 
dt and dt, and dielectric constants e, and e^, respectively (see Fig. 7). 

With the aid of Gauss law, the potential energy difference between layers, A = e{(pb — <Pt), is readily expressed via the gate 
potentials V/ and V/,: 


A = ed 


V, 1 / _U 

d dt / £t db / £b 


(15) 


where Eq is the vacuum permittivity. Assuming that far from gates graphene is undoped (/r = 0), and the source contact is 
grounded, we obtain an expression for the chemical potential in the source region: 


2 d dt / £t db / £b 


(16) 
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Figure 7. (A) Illustration of the electrostatic calculations. Electric field produced by top and bottom gates with potentials V, 
and Vb induces potentials (pi and (pi, at top and bottom layers of GBL. Panel (B) shows the band diagram in the source region. 
Near the grounded source contact there is no bandgap and GBL has zero potential, while under the ’’doping” gate GBL 
exhibits a bandgap A and has potential (tp, + <Pb) /2, which leads to a corresponding midgap shift 


Equations (8)-(16) form a system of equations used for a self-consistent evaluation of A and /r in the source region. The 
system, which can be written in the symbolic form 

(^^)=F(A,m), (17) 

was solved numerically using an iterative method similar to the nonlinear successive over relaxation:^ 

( ^ , (18) 

where i is the iteration number, and X is the relaxation parameter adjusted to achieve convergence. The same approach was 
applied for the evaluation of A and /r in the drain region (with an obvious shift of electric potentials by the drain voltage). 
Considering the GBL band structure under the control (middle) gate, we assumed that the conduction band electrons are in 
equilibrium with the drain contact, while the holes are in equilibrium with the source contact. This approach is justihed as far 
as the tunneling probability is small and the electrons injected from the valence band do not distort the potential distribution. 
To this end, two different chemical potentials for electrons and holes were used in Eqs. (13), (14) when calculating the charge 
density under the control gate. 


2 Distribution of eiectric potentiai, barrier transparency, and current 

The calculation of the tunnel current requires the knowledge of the barrier transparency which, in turn, depends on the distri¬ 
bution of electric potential in the tunneling region. Various studies of the tunnel diodes confirm that it is the maximum field 
in the junction region that governs the barrier transparency.Thus, for an estimate of the tunnel current, the full knowledge 
of potential distribution is not necessary. 

To obtain the junction held, we have solved the Laplace equation in the region between two gate electrodes: the ’’doping 
gate” above the source (held at potential Us) and the top control gate above the channel (held at potential Vg) treating them as 
two keen metal plates.® The lateral distance between the gates is denoted by dg. The effect of the bottom gate on the junction 
held is weak provided db^ dt. The difference between dielectric constants of top and bottom dielectric layers, e, ^ Eb, can 
be approximately taken into account by replacing the top gate oxide thickness with effective oxide thickness, dg/f = dtSb/St- 
The distribution of electric potential in the geometry considered is obtained with conformal mapping technique, and the result 
is 


(p{x,y) 


Us + Vg 
2 


2 
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Vg-Us 

2 


Re 


arcsin 


x + iy'\ 


(19) 


(see Eig. 8), where the origin of the coordinate grid is located between the two gates and lies in their plane. 
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Figure 8. Distribution of the electric potential produced by two keen planes held at potentials Us and Vg- Potential is 
encoded by colour, increasing from blue to red, while arrows represent field lines 


The maximum electric field in the graphene plane (y = —deff) is achieved at x = 0 and equals; 
Vg-Us 


F = 




( 20 ) 


Having obtained the junction field F, we calculate the WKB transparency of the potential barrier ^ {p^,E) separating the 
conduction and valence bands. Assuming that the potential energy depends on x - coordinate linearly, we find 


^{P±,E) =exp 


J Imp|j (x) dx 


XI 


= exp < 


^2 _ ’ 

J Im ^[E + eFx, A (x)] — p^dx 


( 21 ) 


Xi 


where p|| and px are the components of electron quasimomentum that are parallel and perpendicular to the direction of 
junction field, respectively, E is the electron energy measured from the midgap in the source region, the integrals are taken 
between the classical turning points xi and X 2 - Generally, the interlayer asymmetry parameters in the source (A^) and in the 
channel (Ac) are different. For simplicity, we perform a linear interpolation between As and Ac in the junction region. 


A (x) = A 5 + 


Ac —As 
-;- X 


As + 


Ac - A5 

- rx. 

(pc-(ps 


( 22 ) 


where (ps and (pc are the electric potentials in the source and channel regions (obtained in the previous section), and / is the 
length of the transition region between them. 

The tunneling current density Jtun can be written as 

Jtun =egY, ^11 (p) ^ (p) [fs (p) - fc (p)] , (23) 

p 


where V|| = dE /<5p|| is the component of electron group velocity parallel to the current direction, and fs and fc are the Fermi- 
Dirac distributions in the source and channel, respectively. It is more convenient to express all quantities in Eq. (23) via the 
transverse electron momentum p± and its total energy E which are conserved during the tunneling; 

Ey -\-Pmax{E) 

Jy.„=2es/lfs(E)-ME)]^ / (24) 

Ec -PiTiax{£) 
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Figure 9. Illustration of interband tunneling in bilayer graphene. Valence band electrons with small transverse momentum 
and energy far from the band edges (1) pass through the two extrema of the dispersion law while moving in the constant 
electric field and may tunnel from each of these extrema. Electrons with large transverse momentum (2) pass through just 
one maximum. However, during their motion in the forbidden region, there are also two tunneling paths. Finally, for the 
energies close to the valence band edge in Eq. (24), , there are two electron states with positive group velocity 

[(3) and (4)] which contribute to the current independently 


where fs (E) = f{E,jj.s), fo {E) = /{E^iid), fls and fio are the chemical potentials in source and drain regions, respectively 
(we assume the electrons in the conduction band of the channel region are in equilibrium with the drain contact), pmax {E)is, the 
maximum transverse momentum allowed for given total electron energy E in the valence band of source and in the conduction 
band of the channel. More precisely, pmax(E) = min{pc(E),Pv(E)}, where Pc{E) and pv{E) are the inverse functions of the 
electron dispersion law E{p) in the conduction band of the channel and the valence band of source, respectively.^ Finally, 
we note the appearance of the factor of two before the whole expression for the tunnel current which is unique to the band 
structure of graphene bilayer. It appears due to the presence of two points with zero group velocity in the dispersion laws for 
electrons and holes, which implies that a carrier incident on a barrier attempts to pass through it twice (see the red arrow in 
Fig. 9). 

Thermionic currents Jg and 7/, carried by electrons and holes, respectively, are calculated in a similar way, except for 
the absence of the factor of two and barrier transparency equal to unity. The integration limits E in the expressions for the 
thermionic current are restricted to the energies above the barriers created by the doping gates: 

+00 

Je = J Pmax (E) [fs (E) - fo (£)] dE, 

E? 

Jh = / Pmax (E) [fs (E) - fD (£)] dE, 

where Eg and E^ are the edges of conduction band in source region and valence band in drain region, respectively. Electrons 
from the ”hat” (i.e. with p < pmin) do not contribute to the thermionic current because they can neither remain in the ”hat” in 


(25) 

(26) 
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Figure 10. Exact and approximate imaginary part of electron longitudinal momentum in the barrier. Electron with given 
energy and transverse momentum exits valence band at point xi, tunnels through the bandgap, and enters conduction band at 
point X 2 , as shown in the inset 

all three regions (source, gate and drain) nor get out of it (it is possible only at band edges). 

3 Derivation of approximate anaiytic expressions for tunnei and thermionic currents 

Eirstly, we need to obtain an analytic expression for barrier transparency. The integrand in Eq. (21) is too complicated to carry 
out the integration exactly. However, using the fact that the integrand has an arch-like shape (it turns to zero at xi and X 2 and 
reaches maximum somewhere between), we can approximate it with a half ellipse of the same width and height (see Eig. 10). 
In the simplest case of constant A = A,„« the maximum of the integrand is reached in the middle of the barrier and equals 



(27) 


The main contribution to the tunnel current is from electrons with p± < pmin, because for them barrier is the lowest and 
the shortest. Such electrons tunnel from the edge of the valence band to the edge of the conduction band, which gives the 
following barrier length: 


Equations (21), (27) and (28) yield the following expression for barrier transparency in the ’’semielliptical” approximation: 



(29) 


which does not depend on electron energy. To be able to perform the integration with respect to px in Eq. (24), we expand the 
exponent in Eq. (29) in a Maclaurin series up to the second order: 



(30) 


where 



(31) 
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Top gate voltage, Vq (mV) 

Figure 11 . Comparison of gate transfer characteristics at different drain bias obtained via numerical calculations (solid 
lines) and via approximate analytic equations (dashed lines). In the proposed GBL TFET tunneling occurs through varying 
bandgap, so we chose A,j,„ = Aj in order to partially compensate the increase in current introduced by integrating up to 
infinite transverse momenta 


Po = 


I 2 V 2 


a 




Tt \/a- 1 Ar, 


V 


and a = 



2 

. To avoid special functions, we extend the integration region up to infinite p± and obtain 


egkT 




- fs{Ev)fD{Ecy 
JsiEc)fD{Ey_ ■ 


(32) 


(33) 


Now we move on to the calculation of the thermionic current. Near the band edges the dispersion law is approximately 
parabolic: 


Pmax(£') ~Pmin + 


Ayi 

2v2pm 


E 1 
^A"2’ 


(34) 


where plus sign corresponds to the conduction band and minus sign corresponds to the valence band. Taking into account that 
in Eq. (25) the Eermi-Dirac distribution almost coincides with the Boltzmann distibution, we obtain 


Je 


Jh 


egkT 

egkT 



[fs{E^,)-fo{E^.)], 

[/,(Ef)-/z,(Ef)], 


(35) 

(36) 


where index ”S” corresponds to source, and ”D” corresponds to drain. 


4 Comparison with TFETs based on 2d semiconductors with paraboiic bands 

In this section, we compare the dependencies of interband tunnel current on the band overlap in GBL and an equivalent 2D 
semiconductor without the ’Mexican hat’ band structure. The equivalence of the materials implies equal values of band gap 
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Energy E (meV) 

Figure 12. Energy dependencies of electron momentum in the overlapping conduction and valence bands in a 
voltage-biased tunnel junction. Red curves are for graphene bilayer (A = 0.3 eV), blue curves are for an equivalent 
parabolic-band 2D semiconductor. The filled area under the curves is proportional to the interband tunnel current. 


and barrier transparency in the electric field, while the band structure outside the gap is different. For the sake of analytical 
traceability, we assume the barrier transparency to be weakly dependent on electron energy and gate voltage. The low- 
energy band structure of GBL is given by Eqs. (8). The electron dispersion of a parabolic-band material is 

Ec{p)=Eg/2 + ^ (37) 

LlTlc 

in the conduction band, 

E,.{p) = -Eg/2-^, (38) 

Inly 

in the valence band. Here me and iriy are the effective masses for in the respective bands. The interband tunneling occurs 
mainly between conduction and light hole bands as the batTier transparency is small for heavy carriers. 

Similarly to the derivation of Eq. (2) of the main text, one can write down the current density in a parabolic-band 2d 
semiconductor 


Jt 


par 


Ey Pmax{E) Ey 

I dE [fy{E)-fe{E)] J 2dp -2e^% J Praa^{E)dE, 

Ec 0 Ec 


(39) 


where pmax{E) is the maximum transverse momentum of electron with energy E in the conduction and valence bands. The 
value of pmax{E) is the minimum of electron momenta at given energy in the conduction and valence bands; 


Pmax (E) = min {pc{E),py{E)} = min | a/2otc (E-Ec), yj2mv (Ev-E)! • 


(40) 


Graphically, the current 7t,par is proportional to the area under the inverse dispersion laws in the conduction and valence 
bands, this region is filled with blue in Fig. 12. The integral in (39) is evaluated analytically with the result 


7t 


Ay/2 


gs 


par ' 


3 


metny 


nir - 


■ ruv 


■&o{Ey-Ecf\ 


(41) 
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In graphene bilayer, the tunneling current density is given by 

Ey 

Jrm, « 2e^ (2^o) jPm^AE)dE, (42) 

Ec 

where an extra factor of gy = 2 accounts for the valley degeneracy, and the factor of two in front of the barrier transparency 
is due to the two points in electron dispersion Ec{p) with zero group velocity where an incident electron attempts to tunnel 
from the valence band. Again, pmax (E) = mm{pc (E) ,py (is)}, where pc,v (E) are the dependencies of electron momentum 
on energy in the conduction and valence bands. Inverting the electron dispersion, one finds 


Pc{E) = 

Pv{E) = 


1 


\ 


(£-&P + ^ +J(£-£J2(lf + A2)- 


•\ 


(£, - £)2 + ^ + J (£. - £)'()f + Al) - 


4 ’ 

(43) 

4 

(44) 


The energy integral (42) for GBL cannot be evaluated analytically. However, the area under the dependence pmax {E) is 
naturally larger for graphene bilayer than for a typical narrow-gap semiconductor, not to mention an extra factor of four due 
to valley degeneracy and twofold tunneling. At small energies close to Ec, Pc{E) ~ pmin > \/2.mc (E — Ec), while at large 
energies Pc{E)k.E jv, which again lies above \/2mc {E — Ec). In Fig. 5 of the main text, we numerically compare the current 
densities for GBL and 2d semiconductor with parabolic bands using Eqs. (42) and (41). At the operating voltage of 150 mV, 
the current in GBL is 15 times larger than the current in a parabolic band semiconductor, of which the factor of 2 is due to 
the extra valley degeneracy, the factor of 2 is due to the two-fold tunneling, and the factor of 3.5 is due to the ’Mexican hat’ 
dispersion. 


5 Estimate of the band tails in graphene bilayer 

The abruptness of tunnel switching is limited by the tunneling from the localized states in the band tails. ^ This mechanism can 
be alternatively interpreted as tunneling from the band tails, produced by the random dopant-induced potential fluctuations. 
In this section, we evaluate the density of states in graphene bilayer due to these fluctuations. In the quasi-classical limit, the 
’true’ density of states p{E) is evaluated by integrating the bare density of states shifted by eV, po(£ — eV), timed by the 
probability of voltage fluctuation of magnitude eV 


In the above equation, (V^) is the root-mean square voltage fluctuation. 


(46) 


where Vb(q) is the Fourier transform of the potential produced by a single impurity and n, is the impurity density. To obtain a 
non-divergent result for fluctuations in two dimension, a finite distance from the charged impurities to the 2d plane has to be 
considered,^ in this case the Fourier transform of potential is 


Vo(q) 


Kqeiq) 


(47) 


Above, K is the background dielectric constant and £{q) is the dielectric function of graphene bilayer itself. For analytical 
traceability, the latter is taken in the Thomas-Fermi screening approximation:'® 


e{q) = l+qs/q, 


(48) 


where the screening wave vector is calculated as 
= PiE)f{E)[\-m]dE. 


(49) 
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Figure 13. A model band diagram of graphene bilayer - dielectric - metal system used for the evaluation of gate leakage 
current. Graphene bilayer is modelled as a delta-well with potential U{z) = —a5{z-\-dt), where a is chosen to provide a 
correct work function from graphene to the insulating material. 


The integral (46) is evaluated analytically with the result 

{V^) = 2%n,(^^\[qA), (50) 

where K (x) = —e^^{2x+ l)Ei(—2x) — 1. Equations (45) and (50) were used to calculate the density of states in the Eig. 5 of 
the main text numerically. The evaluation of DoS is, actually, an iterative procedure as the screening wave vector depends on 
the DoS itself. However, it converges quickly and the first approximation already is close to the final result. 


6 Modeling of the gate leakage 

To estimate the leakage current from graphene bilayer to the metal gate, one has to multiply the charge density in graphene 
by the area of the gate and by the characteristic ’tunneling frequency’ V?. The latter can be derived from a lifetime of electron 
state bound to graphene bilayer in the presence of continuum of states representing the metal gate (Eig. 13). 

The highest voltages are applied to the ’doping’ gates, and the largest leakage is expected therein. An example of the 
band diagram corresponding to source doping gate is shown in Eig. 13. Graphene bilayer is modelled as a delta-well U{z) = 
— a5{z + dt). The parameter a = 2h y/Ub/2m* is chosen to provide a correct value of the work function from graphene bilayer 
to the surrounding insulator material, m* is the electron effective mass in the insulator. 

A simple wave function matching procedure for the potential profile of Eig. 13 leads to the dispersion equation for the 
quasi-bound states in graphene bilayer 

{k - Kb) {k - ikf) = Kb {k + ikp) , (51) 


where Kb = ^/2mUb/h is the decay constant of the wave function, kf is the electron wave vector in metal (approximately equal 
to the Eermi wave vector), and K = y/—2mE jh is related to the sought-for state energy E <0. The solutions of Eq. (51) are 
complex, E za —Ut, + ihVt/2, where 


V, = 


8t/i, Kskp 

h K^+kjr 


exp[—2 k'b^/] . 


(52) 


Eor non-rectangular barriers, the exponent e should be replaced with the quasi-classical transparency of the barrier 
separating graphene bilayer and the gate. The full leakage current density (per unit width of the channel) is evaluated as 


Jg = eLgV,{ne+nh), 


(53) 
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where «/, is the density of holes induced at the source doping gate (4.5 x lO'^ cm^^ for 11$ = —0.6 V) and rie is the electron 
density at the drain doping gate (8 x lO'^ cm^^ for Ud — 0.25 V). To provide minimum tunneling frequency, we choose 
zirconium oxide as a gate dielectric due to its relatively large effective (tunneling) mass m* = Q3mo and small electron affinity 
XzrOj = 1 -6 eV leading to a quite large value of Ub = XCr - ZzrOz = 2.9 eV. 
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